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Abstract 

The structure of social relations is fundamental for the construction of plausible 
simulation scenarios. It shapes the "way actors interact and create their identity 
■within overlapping social contexts. Each actor interacts in multiple contexts 
within different types of social relations that constitute their social space. In this 
article, we present an approach to model structured agent societies with multiple 
coexisting social networks. We study the notion of context permeability, using a 
game in which agents try to achieve global consensus. We design and analyse two 
different models of permeability. In the first model, agents interact concurrently 
in multiple social networks. In the second, we introduce a context .switching 
mechanism which adds a dynamic temporal component to agent interaction in 
the model. Agents switch between the different networks spending more or 
less time in each one. We compare these models and analyse the influence of 
different social networks regarding the speed of convergence to consensus. We 
conduct a series of experiments that show the impact of different conhgurations 
for coexisting social networks. This approach unveils both the limitations of 
the current modelling approaches and possible research directions for complex 
social space simulations. 
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1. Introduction 


In most simulations of complex social phenomena, agents are considered to in¬ 
habit a space in which structure is very simple. This space has little resemblance 
with the social world it is supposed to depict and for which conclusions are sup¬ 
posed to be extrapolated. This simplicity does not come by chance, rather, it 
is necessary and desired by the researchers: the problems to be approached are 
themselves so complex that whichever factors of complexity can be reduced (or 
at least postponed), the reduction is always welcome. So, it is common practice 
that geographical space is reduced to a two-dimensional grid, and all social re¬ 
lations between agents are condensed into one more or less structured abstract 
relation. Most social simulation and modelling approaches disregard the fact 
that we engage in multiple social relations. Moreover, each kind of social rela¬ 
tion can possess distinctive characteristics that include: rich information such 
as degree of connectivity, centrality, trust, interactions frequency, asymmetry, 
and so on. 

To explore the addition of multiple relations and its consequences for the 
dissemination of phenomena in social simulations, we put forward to reduce 
the emphasis given on agent individual interactions. We accomplished this by 
choosing a simple - and especially neutral - game to model those interactions. 
Our main concern was that the game itself would not provide biases or trends 
in the collective phenomenon being studied, so we chose the consensus forma¬ 
tion game, with a straightforward majority rule to decide the outcome of each 
individual move. 

While most models are simplified descriptions - reductions - of real-world 
phenomena, many constitute complex systems themselves and thus, we need 
techniques such as simulation to explore their properties. We can make a 
distinction between two types of complexity: the complexity of real systems 
{ontological complexity) and the complexity of models {descriptive complex¬ 
ity) [T]. The levels of abstraction in a simulation model can then range from 
data-driven paradigms to more abstract descriptions that allow us to create 
what-if-scenarios. The study of this abstract descriptive complexity in simu¬ 
lation models is as valuable as its data-driven counterpart. We may not be 
able to make predictions based on the direct application of such models; but, 
their study can inspire the engineering of artihcial complex systems and reveal 
properties with applications beyond the explanation of observable phenomena. 
Studying consensus formation can thus advance our understanding of related 
real-world social phenomena. 

Examples of thematic real-world phenomena in social simulation models in¬ 
clude for instance the joint assessments of policies or, in the context of economics 
and politics, the voting problem. Herbert Simon’s investigations on this problem 
were one of the first stepping stones to the field of social simulation [2]; which 
then inspired models (like the one we present in this article) related to evolution 
and dissemination of opinions, which we call opinion dynamics. In agent-based 
opinion dynamics, agent interactions are guided by social space abstractions. 
In some models, the dimensionality and structure of this space is irrelevant 
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(any agent can interact with any other agent). Other models use an underlying 
artefact that structures agent neighbourhoods. Axelrod for instance, represents 
agent neighbourhoods as a bi-dimensional grid in its model of dissemination 
of culture [5]. In an attempt to mimic real social systems, one can also make 
use of complex network models to create the infrastructure that guides agent 
interaction (see [1] for an example). 

In real-world scenarios, actors engage in a multitude of social relations dif¬ 
ferent in kind and quality. Most simulation models don’t explore social space 
designs that take into account the differentiation between coexisting social net¬ 
works. Modelling multiple coexisting relations was an idea pioneered by Peter 
Albin in but without further development. The process of interacting in these 
different complex social dimensions can be seen as the basis for the formation 
of our social identity [ill]. 

In this article, we explore the of modelling opinion dynamics with multiple 
social networks. We look at how the properties of different network models in¬ 
fluence the convergence to opinion consensus. We present a series of models that 
use these networks in different ways: which create distinct emergent dynamics. 
We want to study the consequences of using multiple social networks at the 
same time while maintaining the interaction model as simple and abstract as 
possible (following the methodology in most opinion dynamics literature). We 
present two models where agents can: interact at the same time in the multiple 
networks (choosing partners from any network); or switch between networks 
(choosing only partners from their current network). 

1.1. Article Structure 

This article is organised as follows. In the next section, we present the work 
related to opinion dynamics, social space modelling and complex network mod¬ 
els. In the following section, we describe our game of consensus and introduce 
our multiple model variations. These are designed to study the notion of con¬ 
text permeability and consensus formation in multiple social networks. Section 
1^ describes the experimental setup; the set of tool and methodologies followed 
to conduct our investigations. In section we present and discuss our results 
and compare the different simulation models. Finally, we summarise what we 
learned from our experiments and point out future research directions. 

2. Related Work 

In this section, we present work related to opinion dynamics, social space mod¬ 
elling, and complex network models. 

2.1. Opinion Dynamics and Consensus Formation 

Formal opinion dynamics models provide an understanding, if not an analysis, of 
opinion formation processes. An early formulation of such models was created 
as a way to understand complex phenomena found empirically in groups [5]. 
The work on consensus building (in the context of decision-making) was first 
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studied by DeGroot [3] and Lehrer m- Empirical studies of opinion formation 
in large populations have methodological limitations. Computational sociology 
arises with a set of tools - simulation models in particular - to cope with such 
limitations. We use multi-agent simulation (MAS) as a methodological frame¬ 
work to study such social phenomena in a larger scale. Most opinion dynamics 
models make use of binary opinion values [HI III], or continuous values mill]. 
For a detailed analysis over some opinion dynamics models, refer to m- 

Opinion dynamics models allow us to discover under which circumstances 
a population of agents reaches consensus or polarisation. Agent-based models 
can have broader application outside social simulation though. In computa¬ 
tional distributed systems, consensus is a means by which processes agree on 
some data value needed during computation. Typically, this agreement is the 
result of a negotiation process (often with the aid of a mediator). In human soci¬ 
eties, social conventions emerge to deal with coordination and subsequently with 
cooperation problems m- These conventions are regularities of behaviour that 
can turn normative if they come to be persistent solutions to recurrent prob¬ 
lems. MAS are also capable of producing emergent conventions in coordination 
or cooperation problems HZ]. The decentralised nature of these computational 
models of consensus is highly desirable for dynamic control problems. In these 
scenarios, creating conventions before hand (off-line), or developing a central 
control mechanism for generating them, can be a difficult and intractable task: 
either due to the uncertainty and complexity associated with the environment, 
or due to the system scale and heterogeneity (which makes it very difficult to 
handcraft each component). 

Consensus models are also investigated as means to create conventions in 
MAS. Shoham and Tennenholtz create a bridge between economic literature 
and machine learning by studying a series of models in which multiple agents 
are engaged in learning a particular convention m- They call this process co¬ 
learning. The complexity of these systems comes from its concurrent nature: 
one agent adapts to the behaviour of another agents it has encountered, these 
other agents update their behaviour in a similar fashion, which results in a 
highly non-linear system dynamics. In their work, Shoham and Tennenholtz 
define the notion of stochastic social games and compare different rules with 
the objective of establishing conventions in a decentralised fashion. In their 
stochastic games, they present one particular rule that we use in our models: 

Definition 1. External Majority (EM) update rule: adopt action i if so far it 
was observed in other agents more often than other action and remain with your 
current action in the case of equality. 

This EM rule was shown to coincide with Highest Cumulative Reward (HCR), 
which is a simple rule that states that “an agent should adopt an action that 
has yielded the highest cumulative reward to date.” We use the EM rule in our 
simulation models for its simplicity and success in the evolution of conventions. 
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2.2. Social Structure in Simulation Models 

The usage of a bi-dimensional grid to represent abstract social spaces in simula¬ 
tion models, is one of the most widely used approaches in the agent-based simu¬ 
lation and modelling literature. This has its origins in the “eheckerboard” model 
introduced in computational social sciences by Schelling and Sakoda 
Another famous example that uses bi-dimensional grids as social structure is 
the social simulation model of dissemination of culture from Axelrod |3]. 

Cellular automata (CA) are an example from the area of artificial life that 
also makes use of bi-dimensional structures to model neighbourhood interac¬ 
tions. While these models are idealised frameworks, their exploration can bring 
us deeper insights than models with high level of descriptive complexity (which 
would render their analysis very difficult). The work of Flache and Hegsel- 
mann relaxes the standard CA assumptions about the regular bi-dimensional 
grids by considering irregular grid structures (Voronoi diagrams) [21]. They 
present results on the robustness of some important general properties in re¬ 
lation to variations in the typical grid structure. There are also contributions 
that incorporate both the complexity of network models and the behaviour of 
dynamic processes. One example is the work in |22] . which explores a graph- 
based cellular automaton to study the relationship between spatial forms of 
urban systems and the robustness of different process dynamics under spatial 
change. Moreover, this shows how we can use real geographical information 
to construct graph-based cellular automata (thus making a connection between 
purely abstract models and data-driven models). 

Finally, another type of social space models are the random graph/network 
models. Each complex network - or class of complex networks - captures spe¬ 
cific topological properties. These properties are found in real-world network 
structures. One of the first complex network models is the random graph from 
Erdos and Renyi [ 53 ] . Other examples include the model of preferential attach¬ 
ment [ 21 ] , and the small-world networks [ 53 ] . One common issue of these models 
is the fact that once generated, the network structure is static. The models are 
not really suitable for the description of highly dynamic groups or communities. 
An example of a model that takes such dynamics into account is the model of 
team formation in [26] . A more in-depth review of network models (and their 
properties) is given in [ 57 ] . 

3. Multi-context Models 

In this section, we present our modelling approach to explore the concept of 
permeability between contexts. We use multiple social networks to represent 
the complex social space in which an agent is inserted. In a simulation model, 
this setting can be seen as a n-dimensional scenario where each dimension con¬ 
tains a network that represents a different social relation (see figure . The 
word “context” is used in many different senses and it is by itself subjected to 
many analysis [28] . Here, we use the term social context in a simpler and more 
strict way. Agents belong to distinct social contexts which are their neighbour¬ 
hoods in these multiple networks. The social context of an agent is thus a set 
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of neighbours currently available for interaction at a particular network or set 
of networks. Context permeability is the ability of a particular norm, strategy, 
opinion, or trend, to permeate from one social relation to another. Agents be¬ 
longing to disjoint neighbourhoods in different networks serve as bridges between 
subsets of the population that wouldn’t otherwise be influenced by a consensus 
formed in a distant coalition. 



Figure 1: Multiple social network structure that shapes the social space 
in our simulation models. 


We study the notion of context permeability in different simulation models 
in which agents interact using a simple consensus game. The society of agents 
has to adopt a binary opinion value according to a majority rule. We use the 
speed of consensus as a measure for self-organisation and explore the relationship 
between different network topologies and this measure. 

We partitioned the design space of our experiments using different simu¬ 
lation models. Each model is designed to analyse different aspects of context 
permeability. In a first model [liisn], we study the notion of permeability by 
overlapping social networks. In a second model |12j . we analyse the dynam¬ 
ics introduced by switching between social contexts (neighbourhoods in different 
networks). This adds a temporal component that changes the dynamics of the 
game: agents can only perform encounters with neighbours that are active in 
their current context. (It also adds the chance for some agents to be isolated 
during the simulation.) 

3.1. Simulation and Consensus Game 

Our agent-based models are designed as discrete-event models. On each simula¬ 
tion cycle, every agent executes a simulation step. The agents are selected in a 
random uniform fashion for execution. This is common practice in agent-based 
simulation and guarantees that there is no bias caused by the order in which 
the agents are executed. For each simulation model, we present a description 
of the individual agent behaviour. Every model we present can be seen as a 
binary opinion dynamics model, or a consensus game. In this game, the soci¬ 
ety of agents tries to reach an arbitrary global consensus about two possible 
choices or opinion values. During a simulation run, each agent keeps track of 
the number of each opinion value of their interaction partners. In each itera¬ 
tion, each agent selects an available neighbour to interact with and observes its 
current opinion value. The agent decides to switch its current opinion choice 
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if the observed opinion value becomes the majority of the two possible choices. 
The goal of the game is to reach a global consensus, but the particular choice 
that gets collectively selected is irrelevant. What is important is that overall 
agreement is achieved. In the consensus game we consider, each agent uses the 
previously presented external majority (EM) rule to update its opinion value 
(see definitionin section]^. 

The reason we are not interested in which exact option gets selected is the 
same reason why we chose such a content-neutral game for the individual inter¬ 
actions. From the very beginning, our research questions were focused on the 
collective properties, the structure of the networks involved, and the dynami¬ 
cal outcomes of the simulations: thus leaving out issues related with individual 
motivations and desired results (either for any individual agent or the society). 
Granted, most applications of this research will often include such goals, for 
inquisition into complex social issues, including possibly policy design, in which 
the individual rationality must be considered, and will have a very strong influ¬ 
ence on the collective outcomes. 

3.2. Context Permeability 

The first model of context permeability Hi ISO] is designed to study the setting 
where an agent is immersed in a complex social world: where it engages in 
a multitude of relations with other agents. In this model, we consider that 
the context permeability is created by agents that can belong simultaneously to 
several relations. As an example, two agents can be simultaneously family and 
co-workers. While some links might not connect two agents directly, others can 
contain neighbourhoods in which they are related (either directly or through a 
common neighbour). Model describes the behaviour of each agent for each 
simulated step. 
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Model 1 Context Permeability 


M {Simulation Model / Environment} 
A (Current Agent in Execution} 


(Randomly select a network} 

r •<— randomUniform(0,M.netu'orfcs 
network ^ M.networks[r\ 

(Randomly select a neighbour} 

neighbours network.neighboursOf(A) 
r •<— randomUniform(0,neighbours.length()) 
partner ^ neighbours [r] 

(Update opinion and memory} 

(opinions take the values {0,l}~so they can be used as an index} 
A. memory [partner. opinion] H—h 

if (A.memory [partner.opinion] > A.memory [A.opinion]) then 
A.opinion ^ partner.opinion 

end if 


In this model, on each time step, each agent select a random neighbour from 
any social network available and updates its opinion based on the opinion of the 
selected partner using the external majority rule. 

3.3. Context Switching 

In the context switching model |12] . agents are embedded in multiple relations 
represented as static social networks and they switch contexts with a probability 
(jCi associated with context Ci. The switching probability is the same for every 
agent in a given network. Each agent switches from its current neighbourhood 
to another network. The agents are only active in one context at a time, and 
can only perform encounters with neighbours in the same context. We can think 
of context switching as a temporary deployment in another place, such as what 
happens with temporary immigration. The network topology is static; when 
an agent switches from one network to another, they become inactive in one 
network and active in their destination. 






Model 2 Context Switching 

M {Simulation Model / Environment} 
A (Current Agent in Execution} 


(Get the current network {context is a network index)} 

cNetwork ^ M.networks[A.context] 

(Get the neighbours for the current network} 

neighbours ^ cNetwork.neighboursOf(A) 

(hlter by agents active in the same network} 

neighbours ^ (n | n € neighbours A (n.context = A.context)} 

(Randomly select a neighbour} 

r ^ randomUniform(0,neighbours.length))) 
partner •<— neighbours [r] 

(Update opinion and memory} 

(opinions take the values {0,l}-so they can be used as an index} 

A. memory [part ner .opinion] H—h 

if (A.memory [partner.opinion] > A.memory [A.opinion]) then 
A.opinion ^ partner.opinion 

end if 

{Switch to another network} 

switchingProb ^ M.params.switchingProb(A.context) 
r ^ randomUniform(0,l) 
if (r < switchingProb) then 

numNets ^ M.networks.length)) 

nextNetworks ^ (i | )i € [0,numNets[) A ) M. net works [i] != A.context)} 
r •<— randomUniform)0,nextNetworks.length))) 

A.context ^ nextNetworks[r] 

end if 


In this model, agents select an )active) partner from their current neighbour¬ 
hood. If a partner is available, the agents update their opinion value based on 
EM. At the end of the interaction, the agent switches from the current context 
to a different one with a probability associated with its current context. 

This model describes an abstract way to represent the time spent on each 
network using the switching probability ^<7.. Here, the permeability between 
contexts is achieved using this temporal component. Context switching intro¬ 
duces a notion that has not been explored in the literature so far: the fact that, 
although some social contexts can be relatively stable, our social peers are not 
always available for interaction and spend different amounts of time in distinct 
social contexts. 
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4. Experimental Setup 

In this section, we present all the tools and processes necessary to produce 
the current research output. The experiments were developed using the MA¬ 
SON [3T] simulation framework written in Java. All the code for the simulation 
models can be found here [52]. Moreover, all the results presented in this article 
were made reproducible by using the statistical computing language R [33j and 
the R package knitr |34j . Knitr is used for generating reports that contain the 
code, the results of its execution (plots and tables), and the textual description 
for such results. All the R code used to produce the data analysis and the 
conhguration files used to reproduce the data can be found in [35j . 

For each conhguration in our experiments, we performed 100 independent 
simulation runs. The results are analysed in terms of value distribution, average 
and standard deviation, or variance over the 100 runs. Like we described in 
section [3T] our models are discrete-event models. In each cycle agents execute 
a step in a random order: this is so that this dynamical system -our simulation 
model- is not inhuenced by the order in which agents are executed. Simulations 
end when total consensus is achieved or 2000 steps have passed. 

^.1. Social Network Models 

The social networks used in our models are k-regular and scale-free. A k-regular 
network is a network where all the nodes have the same number of connections. 
These networks are constructed by arranging the nodes in a ring and connecting 
each node to their next k neighbours. Each network has 2k edges per vertex. 

Scale-free networks are networks in which the degree distribution follows 
a power law, at least asymptotically. That is, the fraction P{k) of nodes in 
the network having k connections to other nodes goes for large values of k 
as: P{k) ^ k~'^ where 7 is a parameter whose value is typically in the range 
2 < 7 < 3. 

We use the method proposed in [21j by Barabasi and Albert to construct the 
scale-free network instances using a preferential attachment. This model builds 
upon the perception of a common property of many large networks: a scale- 
free-power-law distribution of node connectivity. This feature was found to be a 
consequence of two generic mechanisms: networks expand continuously by the 
addition of new vertices, and new vertices attach preferentially to sites that are 
already well connected (more commonly known as “preferential attachment”). 
In these networks, the probability P{k) of two nodes being connected to each 
other decays as a power law, following P{k) ^ k~^. 

The network instances were generated using the h-have network library |36] . 
This is a Java library that allows the creation and manipulation of network/graph 
data structures. It also includes the more commonly used random network mod¬ 
els from the literature. 

4 . 2 . Measuring Network Properties 

We measure two structural properties of our networks: the average path length 
and the clustering coefficient. The average path length measures the typical sep¬ 
aration between two vertices in the graph. The clustering coefficient measures 
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the average cliquishness of the graph neighbourhood. This measure quantihes 
how close the neighbours of a node are to forming a clique. Duncan J. Watts and 
Steven Strogatz [25] introduced the measure to characterise a class of complex 
networks called small-world networks. 

These properties characterise structures that can be found in the real scenar¬ 
ios. The small-world networks generated by the models of Watts and Strogatz 
for instance, model real-world social networks with short average path length, 
local clustering structures, and triadic closures between the nodes [25]. One 
shortcoming of this model, is the fact that these networks do not possess node 
hubs found in many real-world social networks. (Nevertheless, these models are 
designed to create networks with specific properties, not as general models for 
all kinds of real social networks) 

Since we construct scenarios where we overlap highly clustered k-regular 
networks and scale-free networks with node hubs, we are interested in the kind 
of properties that emerge from merging these structures - and if/how they are 
comparable to the existing models -. The networks were produced by the b- 
have network library |36j . exported to files, and analysed using the igraph R 
package m 

5. Results and Discussion 

In this section, we analyse and discuss the results for different experiments 
with both the context permeability and context switching models. We start by 
analysing the structural properties that arise from combining different networks. 
We then observe the impact of different model configurations in the convergence 
to consensus. We correlate the outcomes of our context permeability model sim¬ 
ulations with the structural properties of the multi-network structures. We also 
explore how the new dynamics introduced by the context switching mechanism 
affects the consensus building process. 

5.1. Overlapping Network Properties 

In this first analysis, we investigate the properties of different network topologies 
used in our models. One of the parameters in some experiments is the number 
of networks in which the agents interact. Adding more networks implies the 
addition of more connections. We show that what is important is not just the 
the number of connections, but the properties of the resulting structure (when 
we merge the networks). We analyse the networks as follows. 

Each network is generated in such a way that the node indexes are ran¬ 
domised. This means that we can have multiple networks with the same topol¬ 
ogy and each node can have different neighbourhoods in different networks. 
Also, neighbourhoods do not necessarily overlap due to the node shuffling. 
Adding more networks to the social space is not the same as creating a sin¬ 
gle network with twice the number of connections. 
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Figure 2: Overlapping of two k-regular networks, each with k=10. 


Table 1: Properties for overlapping of two k-regular networks, each with 
k=10. 



Nodes 

Edges 

Clustering Coef. 

Avg. Path Length 

10 -regular 1 

100 

1,000 

0.711 

2.980 

10 -regular 2 

100 

1,000 

0.711 

2.980 

Combined 

100 

1,790 

0.512 

1.638 

20 -regular 

100 

2,000 

0.731 

1.788 


Consider figure]^ We created two random k-regular networks with k=10 and 
merged them ignoring edges from common neighbours. (In the simulation mod¬ 
els this edge agglutination is not performed, we use it here for the purposes of 
network analysis.) Common edges are highlighted with a different colour. 

The networks are drawn using the 
Kamada-Kawai Layout [38] which treats 
the “geometric” (Euclidean) distance be¬ 
tween two vertices in the drawing as the 
“graph theoretic” distance between them 
in the corresponding graph. We can 
see by the results of the layout in fig¬ 
ure that the distance between nodes 
in the network has decreased, this is also 
confirmed by the reduced average path 
length of the combined networks (see ta¬ 
ble!^. To illustrate our point, we anal- Figure 3. Single k-regular 
ysed one k-regular network with k = 20 network with k=20. 

(see figure and the corresponding entry on table . The clustering coefficient 
of a 20-regular network is approximately the same as the one of a 10-regular 
network (higher than the combination of two 10-regular networks). The average 
path length also decreases as there are more connections between previously dis¬ 
tant nodes. As we will show, these properties may vary due to the node indexes 
being subjected to random permutations. 
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Since we model the agent social space as a multitude of networks, we need 
to investigate the properties resulting from the merging of these networks. We 
do this empirically by taking the network instances used in the simulations, and 
analysing the distribution of clustering coefficient and average path length for 
the different configurations. 

5.1.1. Properties of Overlapping K-Regular Networks 

We will now investigate what kind of properties we can get from merging k- 
regular networks. First, we analysed the distribution for the average path length 
and the clustering coefficient values over 100 network instances (each with 100 
nodes). From a box plot preliminary analysis (see figures |A.2| and 
see that the average path length does not vary much for the 100 instances, as 
such, the average makes a good descriptor for these properties. 

Figurej^shows the average value for the average path length of 100 k-regular 
network instances (each with 100 nodes). We can see that the average path 
length changes more drastically when we go from 1 to 2 networks. This is 
precisely the effect we can see in figure Merging these networks at random 
effectively creates multiple shortcuts between points that were not connected in 
the original k-regular topology. Beyond this point, adding more networks does 
not modify this structural property in a significant way. Note that, since each 
network has 100 nodes, with fc = 50 the network is fully connected (hence the 
average path length being 1). 


A.l I, we can 




Figure 4: Average value for the average path length of 100 instances 
of overlapping k-regular networks (containing 100 nodes each) with 
k={l,2,3,4,5} pal), and k= {10,20,30,40,50} pbl). 


Figure shows the results for the average clustering coefficient of 100 instances 
of overlapping k-regular networks. Adding multiple networks changes the clus¬ 
tering coefficient of the resulting structure. For k < 20, the clustering coefficient 
drops when we add a second network. This happens because when we merge 
these networks, their connectivity is not high enough to generate shortcuts with 
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tightly clustered neighbourhoods. The clustering coefficient is computed by 
taking the average local clustering coefficient of each node in the network. An 
initial highly clustered k-regular network (with all the nodes having the same 
structure) is affected when you modify the node neighbourhoods in such a way 
that some are more clustered than others. Networks with fc = 1 have a ring-like 
structure, so there are no nodes forming a cliqu^with their neighbours. 



Figure 5: Average value for the clustering coefficient of 100 instances 
of overlapping k-regular networks (containing 100 nodes each) with 
k={l,2,3,4,5} (l5a]), and k = {10,20,30,40,50} (ISbl). 


5.1.2. Properties of Overlapping Scale-free Networks 

We also looked at the properties that result from merging scale-free networks 
(see figure A.3). Like the previous results, the properties don’t vary much 
between the 100 different instances. We can then say that specific configurations 
display a specific average path length and clustering coefficient both for k-regular 
networks and scale-free networks. We plotted the average of the average path 
length (figure 6a I and clustering coefficient (figure A.3b I for each overlapping 
network configuration. 


clique is a group of nodes such that every two nodes are connected by an edge. 
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2 3 4 5 

Number of Networks 


(a) 



Figure 6: Average value for the average path length (A.3a I and clustering 
coejflcient (A.3b) for 100 instances of overlapping scale-free networks 
(containing 100 nodes each). 


The parameter d (figure ^ dictates how many edges are added by preferential 
attachment (see sectior {4T I. For d = 1, the resulting network is a forest: a 
network composed of disjoint tree graphs. Figure shows us that the average 
path length still decreases consistently when we add more networks. One of the 
characteristics of scale-free networks is its short average path length, thus, for 
d > 2, adding more edges does not make a significant difference. Scale-free 
networks also have a reduced clustering coefficient (see figure 6b). When we 
add more networks the clustering coefficient increases. 


5.1.3. Merging K-Regular with Scale-Free Networks 

Finally, we look at what happens when we merge k-regular networks with scale- 
free networks. We didn’t make an exhaustive analysis for this configuration, 
but we show what is the resulting structure of such merging. You can see in 
figure that this results in a network with a lower average path length due to 
the shortcuts created by the scale-free network. It is also less dense than the 
structure with two 10-regular networks (figure]^. We consider only the merging 
of two networks. 
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Figure 7: Overlapping of a 10-regular with a scale-free network with 
D=2. 



(a) 


(b) 


Figure 8: Average average path length for 100 instances of heterogeneous 
network conhgurations: one scale-free with d = {1, 2,3,4, 5} and one k- 
regular network with k = {1,2,3,4, 5} (8a) and k = {10,20,30,40,50} 
(8b). (See the respective box plots in the appendix: hgure A.4). 


Figure shows the average path length of 100 independently generated in¬ 
stances of combinations of one k-regular network and one scale-free network 
with different values k and d respectively. For configurations where the k-regular 
networks have low values of k (figure [8^, adding a scale-free network always de¬ 
creases the average path length. Moreover, the higher the d, the lower the path 
length (as more connections are formed between different nodes). With higher 
values of k (figure |8b| ) however, increasing the value of d does not decrease the 
average path length. This is not surprising given the level of connectivity of this 
structure. The first network is already highly dense and connected, so adding 
a few extra connections doesn’t matter; what does matter is that we have a 
scale-free network creating shortcuts between nodes that were not previously 
connected. 
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Figure 9: Average clustering coefficient for 100 instances of hetero¬ 
geneous network configurations: one scale-free with d = {1,2,3,4, 5} 
and one k-regular network with k = {1,2,3,4,5} (A.3b) and k = 
{10,20,30,40,50} (See the respective box plots in the appending: fig- 
urdA.5|). 


Figure shows the results for the clustering coefficient. These are similar to 
the results we get for k-regular networks (see figure . The difference is that 
we are not adding more networks but rather increasing the parameter d of the 
scale-free (number of connections added by preferential attachment for each 
note). In this case, the effects are qualitatively similar with the ones of merged 
k-regular network with low value of k. In contrast, for high values of k (figure 
9b), a higher value of d is not enough to make the network more clustered. 


This happens because scale-free networks are more sparse and have a very low 
number of triangles between nodes. 


5.2. Context Permeability 

In this section, we discuss the results for various sets of experiments on the eon- 
text permeability model |291130j . First, we used k-regular and scale-free networks 
with our model. Each network layer was configured with the same topology. We 
analyse the context permeability model in terms of convergence to consensus: 
the ratio of simulations that converged to total consensus over 3000 runs; and 
how many agent encounters were needed for this to happen. 


5.2.1. Convergence Ratio 

To analyse the ratio of convergence to consensus for different network configu¬ 
rations, we correlate this ratio with the properties that the combined network 
structure exhibits. As we have seen in section |5.H the properties don’t vary 
much for different random instances, so we use the average as the descriptor for 
these properties. 

The following tables show the results of our simulations. Table shows 
the percentage of total consensus achieved with 2000 cycles for 3000 runs. For 


17 













k-regular networks with a small k, consensus is rarely achieved with a single 
network. However, as soon as we add more networks, consensus is achieved in 
a significantly greater number of occasions. These results are especially inter¬ 
esting for the scale-free networks (table [^: as soon as we add more networks, 
for a scale-free with d > 2, we can achieve consensus in a lot of occasions. This 
is to show that achieving consensus is not just a matter of connectivity. With 
scale-free networks, we can achieve consensus more often with less clustered 
networks. This convergence is achieved not because of the clustering coefficient 
but rather due to lower average path length. This is difficult to visualise us¬ 
ing only the tables and the previous plots so, to confirm this hypothesis, we 
measured the Spearman’s corre/ofzor^ between the ratio of convergences to con¬ 
sensus, the average average path length, and the average clustering coefficient 
for the corresponding network configurations (see table . 

Table 2: Ratio of convergence to total consensus in 3000 independent 
runs with with a number of concomitant networks {ff nets.) equal in 
kind: k-regular networks with A: = {1, 2,3,4, 5,10, 20, 30,40, 50}. 


# nets. 

1 

2 

3 

4 

value for k 

5 10 

20 

30 

40 

50 

1 

0.00 

0.00 

0.00 

0.01 

0.02 

0.18 

0.67 

0.97 

1.00 

1.00 

2 

0.08 

0.59 

0.81 

0.90 

0.93 

0.98 

1.00 

1.00 

1.00 

1.00 

3 

0.66 

0.94 

0.98 

0.99 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

4 

0.88 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5 

0.95 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


Table 3: Ratio of convergence to total consensus in 3000 independent 
runs with a number of concomitant networks {ff nets.) equal in kind: 
scale-free networks with d = {1, 2, 3,4, 5}. 


# nets. 

1 

value for d 

2 3 4 

5 

1 

0.00 

0.29 

0.86 

0.96 

0.99 

2 

0.28 

0.93 

0.99 

1.00 

1.00 

3 

0.78 

0.99 

1.00 

1.00 

1.00 

4 

0.93 

1.00 

1.00 

1.00 

1.00 

5 

0.97 

1.00 

1.00 

1.00 

1.00 


^The Spearman’s rank correlation coefficient is a non-parametric measure of statistical 
dependence between two variables. It assesses how well the relationship between two variables 
can be described using a monotonic function. 
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Table 4: Spearman’s correlation between convergence ratio (CR), average aver¬ 
age path length (APL) and clustering coefficient (CC) along with the respective 
p-value and confidence interval, for a confidence level of 95%. 


network model 

X 

Y 

Correlation 

CI(95%) 

p-value 

k-regular 

CR 

APL 

-0.726 

[-0.836, -0.561] 

2.44 X 10-9 

k-regular 

CR 

CC 

0.264 

[-0.015, 0.505] 

6.40 X 10-2 

scale-free 

CR 

APL 

-0.912 

[-0.961, -0.808] 

2.28 X 10-1° 

scale-free 

CR 

CC 

0.633 

[0.318, 0.823] 

6.78 X 10-4 


5.2.2. Number of Encounters to Achieve Consensus 

We now analyse the convergence in terms of average number of meetings during 
a simulation run. Since the maximum number of simulation cycles is 2000, the 
maximum number of encounters is 2 000 000 (we have 100 agents and each 
one performs one encounter per cycle). We show that for some configurations, 
the average of encounters is not a good descriptor since it varies greatly from 
run to run. Rather than considering the average, it is better to look at the 
distributions of these measures. Nevertheless, we included the data with all the 
average number of encounters, as well as the respective standard deviations for 
different network configurations. You can find this data in appendix [B} 

First, we present the results for k-regular networks. To observe the distri¬ 
bution of the number of encounters, consider the box plots for A: = {1, 2,3,4, 5} 
(figure [TO^ and k = {10, 20, 30,40, 50} (figure [l0b| . Figure [I^ shows the results 
for configurations with a number of networks > 3 (since less networks lead to 
much worse results in terms of convergence speed). These plots reveal two types 
of situations: the typical runs, in which the number of encounters doesn’t vary 
much, and some “outlier runs,” which present some interesting dynamics, as we 
will see later on. 
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Figure 10: Number of meetings over 100 independent simulation runs 


for multiple k-regular networks: k = (1,2,3,4,5} (10a) and k = 


{10,20,30,40,50} (10bI). 
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Figure 10 shows us results qualitatively similar to what we previously observed 
for convergence ratios (see section 5.2.1). Networks with lesser connectivity and 
a greater average path length lead to total convergence less often. The simula¬ 
tion outcomes also vary much more in comparison with the other configurations 
(fc > 1 and networks > 3). 

Notice that for fc = 3 and networks = 3 (10b I, we still have runs that deviate 
from what we consider the typical behaviour. Instead of disregarding this as 
an outlier, we explore this behaviour to better understand the dynamics of our 
External Majority consensus game. 

We analysed the network structure from this run and observed that there 
were no distinct differences between the network instances from the “outlier 
run” and the ones of the “typical” run. This pointed to the consensus game 
itself as a possible cause for the irregular behaviour - not to the network struc¬ 
ture -. What happens is that the agents reach a state from which converging 
towards consensus is significantly harder. Villatoro [51] found that some com¬ 
plex network models lead consensus games to form metastable sub-conventions. 
These are very difficult to break and reaching 100% agreement is not as straight¬ 
forward as assumed by previous researchers. We are not concerned with this 
difficulty, instead, we want to explore how to describe what happens in simi¬ 
lar cases. Figure El shows the consensus progression for two simulation runs a 
typical case and the outlier. 



(a) Typical Run 


(b) Outlier Run 


Figure 11: Number of agents with each opinion during a simulation run: 
typical run in 


11 a 


and “outlier run” in 11b 


We take the outlier in the configuration with fc = 3 and networks = 3, and 
look at the memory that agents use in the consensus game. As we described 
previously, agents record the number of individuals encountered with each opin¬ 
ion value (two possible values in this case). We look at the average difference 
between the number of values observed for all the agents. We also look at the 
variance to find how the opinion observations evolve throughout the simulation. 
Figure [l^ shows the average difference between opinion memory for the first 10 
steps (which was when the simulation converged for the normal run). 
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Simulation Step Simulation Step 


(a) Typical Run 




(b) Outlier Run 


Figure 12: Average and variance of the difference between opinions 
observed for the 100 agents: typical run in figure [T2a] and first steps of 
the “outlier run” in figure 


12 b 


We can see that both the average difference in the opinion memory and 
the variance grow faster in the normal run than in the outlier run. This is 
not enough to establish a significant difference between the two. Moreover, 
the memory differences are quite similar, but the number of agents for each 
opinion was quite even. This happens due to the position of agents in the 
network and their initial opinion values. These circumstances led the agents to 
a initial stability. Note that the opinion difference continues to rise throughout 
the simulation along with the variance. The variance starts to increase more 
rapidly after 200 steps. The exponential growth reveals that the convergence 
to consensus is not done evenly throughout the network, this is one of the 
reasons why consensus takes more time to be achieved. Also, the low variance 
in the first 200 steps shows us that the opinion strength was evenly distributed: 
we have almost the same number of agents for each opinion but the memory 
differences are qualitatively the same for both opinions. As soon as this variance 
increases, one of the opinions breaks the stable state and convergence towards 
global consensus begins. 
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Simulation Step Simulation Step 

Figure 13: Average and variance of the difference between opinions 
observed for the 100 agents throughout the outlier run. 


In other words, different neighbourhoods in the network lead to different opinion 
observation configurations (possibly even creating self-reinforcing structures). 
Thus, the influence of heterogeneous connectivity is reflected on the growth of 
opinion difference variance. 
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Figure 14: Number of meetings over 100 independent simulation runs 
for multiple scale-free networks: with d={l,2,3,4,5} (14a) with a zoom 


on the best configurations in 14b 


Figure [I^ shows the average number of encounters during 100 independent sim¬ 
ulation runs for scale-free networks. (We don’t include the configuration for one 
network and d = 1 because these never converge.) The convergence ratio is very 
low for the configurations with one network d = 2 and 2 networks with d = 1; 
hence the high number of encounters (see also table [^. 

These results are similar to those of k-regular in that adding more networks 
speeds up the convergence to consensus. The difference is that more “atypi¬ 
cal” runs occur with multiple scale-free networks. Again, this is the effect of 
these types of topologies. As we discussed previously, in these less connected 
topologies, agents are more prone to be arranged in such a way that progres- 
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sion towards consensus is more difficult to attain due to bottlenecks. Different 
regions in the network might converge towards different opinion values. 

5.2.3. Heterogeneous Network Configuration 

We also performed experiments with mixed network topologies: we used one 
scale-free and one k-regular networks. Table shows the convergence ratio for 
the heterogeneous network configuration with the multiple values for d and k 
for the scale-free and k-regular respectively. 

Table 5: Ratio of convergence to total consensus in 3000 independent 
runs with two networks: one k-regular and one scale-free network with 
different k and d values respectively. 


d 

1 

2 

3 

4 

k 

5 

10 

20 

30 

40 

50 

1 

0.18 

0.47 

0.55 

0.67 

0.61 

0.76 

0.77 

0.81 

0.75 

0.82 

2 

0.63 

0.72 

0.88 

0.89 

0.90 

0.97 

1.00 

1.00 

0.99 

1.00 

3 

0.77 

0.86 

0.87 

0.97 

0.98 

0.99 

0.98 

0.99 

1.00 

1.00 

4 

0.73 

0.92 

0.92 

0.97 

0.96 

0.96 

0.99 

1.00 

1.00 

1.00 

5 

0.87 

0.85 

0.94 

0.94 

0.98 

0.99 

1.00 

1.00 

1.00 

1.00 


One of the differences between this configuration and the previous ones can be 
observed in the simulations with d = 1 and k = 1. The convergence ratio is 
higher than with 2 1 — regular networks (see table - but not better than two 
scale-free networks. This happens because with k = 1 the network is basically a 
ring and has the maximum possible average path length for a connected graph. 
Adding two rings improves the convergence ratio but the underlying structure 
is still very susceptible to self-reinforcing structures (a connected sub-graph is 
basically a line). Adding a scale-free changes this drastically as we are mixing 
tree-like network components with a ring. 
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(a) 


(b) 


Figure 15: Number of meetings over 100 independent simulation runs 
for heterogeneous configuration: one k-regular network with k = 


{1, 2, 3,4, 5} and one scale-free network with djl, 2,3,4, 5}. Figure 15b 
shows a zoom on the best configurations. 


Figure [TS] shows the results for the number of encounters necessary for the 
agents to achieve consensus. We focus on the k-regular networks with higher 
values of k. We can see that with these two networks, the configurations with 
values of d >= 2 for scale-free networks produce drastically better results in 
terms of speed of convergence. The data for the average number of encounters, 
as well as the respective standard deviations referent to the results in figure [T^ 
can be found in appendixtable |B.3[ 


5.3. Context Switching 

In this section, we present the results relative to the context switching model. 
The difference between this model and the previous one is that the agents no 
longer interact in multiple networks at the same time (being able to select any 
neighbour at a given simulation step). In this model, they become active in a 
single network at a time. Agents can switch to a different network at the end 
of each step. The switching mechanism uses a probability associated with each 
network. With this new idea of swapping contexts, we covered some space left 
undeveloped in Antunes and colleagues original work [30] . 

As we discussed in section 3.3 (see model |^, the switching probability dic¬ 
tates the frequency with which an agent switches from the current network 
after an encounter has been performed. In an abstract manner, this allows us 
to model how much time agents spend on each network. A further development 
of this model will be to assign different preferences to different networks for 
each agent. This can help us model phenomena such as real-world agents that 
dedicate more or less time consuming content from different social networks. 
For now, we attribute this probability to the network to restrain the model 
complexity. 
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5.3.1. Exploring the Switching Probability 

We first look at the influence of the new parameter (switching probability) in the 
convergence to consensus. We span the switching probability between 0 and 1 in 
intervals of 0.05. We do this for 2 networks with k-regular and scale-free topolo¬ 
gies. Figure [T6| shows the results for two k-regular networks: one configuration 
with A: = 10 and one configuration with k = 30. We use these two configura¬ 
tions to study the relationship between the number of encounters, the switching 
probability parameter, and the connectivity of the networks. Remember that 
with fc = 30 the average path length is also lower (section 5.1.1). 



(a) 10-regular Networks 



0.00 0.25 0.50 0.75 


5000 



Switching for Network 1 
regular network, k=30 


Average Encounters 

100000 

75000 

50000 

25000 


(b) 30-regular Networks 


Figure 16: Contour plot for the average number of meets during a sim¬ 
ulation for 100 independent runs: 2 10-regular networks ( |16a[ ) and 30- 
regular networks (16b). See the perspective plots in appendix [C[ 


For smaller values of fc, symmetry in the context switching probability (having 
the same probability in both networks) is more important if the agent switches 
less from one of the networks. Switching less from one network means spending 
more time in that network. This means that switching more from the other 
network can be disruptive to a neighbourhood that has already converged to a 
sub-convention. This can be observed in configurations with k-regular networks 
(figure [16^ but it is especially apparent in scale-free networks (figure 17aI. 
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Switching for Network 1 
scaie-free network, d=1 


(a) Scale-Free d — 1 



(b) Scale-Free d = 5 


Figure 17: Contour plot for the average number of meets during a simu¬ 
lation for 100 independent runs: 2 scale-free networks with d = 1 (17a) 
and d = 5 (17b). See the perspective plots in appendix 


When we increase the connectivity (and lower the average path length), 
both in k-regular and scale-free networks, the switching probability becomes 
“irrelevant” -for any probability value, the speed of convergence will be the 
same-. The exception is obviously when considering the switching probability 
with value 0. In this case, agents are isolated (they are initially equally dis¬ 
tributed throughout the networks) and reaching consensus becomes very hard. 
(Especially because agents are active only in one network and this can create 
disconnected graphs.) 

One major difference from the previous model is that the agents no longer 
interact with a neighbour from any network at any time. The neighbour selec¬ 
tion has to be made from the network they are currently active in. This causes 
the agents to spend different amounts of time in different networks. An agent 
can be very influential in one network - being a very central node in that context 
- and marginally important in another: in a sense its opinion in one network 
can be more important to overall convergence than in other network, due to its 
position in the network. 

To explore the contribution of connectivity to the new context switching 
dynamics, we setup two scale-free networks; one with d = 1 and the other with 
d = 5. We then spanned the switching probability in increments of 0.05. The 
results can be seen in figure 
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Switching for Network 1 
scale-free network, d=5 


Figure 18: Contour plot for the average number of meets during a sim¬ 
ulation for 100 independent runs: 2 scale-free networks, the first with 
d = 5 and the second with d = 1. See the perspective plots in appendix 
[C|(figure[c7). 


In this case it is more important (to convergence speed) to switch from the 
network with lower average path length and a forest-like composition (and a 
probability of at least the same value or more than the one of the other network 
with d = 5). The same can be observed when we mix k-regular with scale-free 
networks (see figure [l^ . 

Switching less from the less connected network is bad in both cases: it is 
possible that spending more time in the network with bigger neighbourhoods 
allows for a stabler convergence. Sub-conventions can emerge in scale-free net¬ 
works with d = 1 because usually some network regions are isolated by a single 
node - the root of the sub-tree they belong to. 
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Switching for Network 1 
reguiar network, k=10 


Figure 19: Contour plot for the average number of meets during a sim¬ 
ulation for 100 independent runs: one 10-regular network (k=10) and a 
scale-free network with d = 1. See the perspective plots in appendix [C] 
(figure C.6). 


5.3.2. Comparison with Context Permeability 

Finally, to compare the context switching model with the previous model of con¬ 
text permeability, we froze the switching probabilities with values C, = 0.25,0.5, 0.75 
and varied the number of networks. (We considered these values because they 
conveniently characterise the probability in low, medium and high switching.) 
We compare the results with the context permeability in terms of average num¬ 
ber of encounters to achieve consensus. Note that in some cases, due to the 
distribution of the number of encounters, the average is not an accurate de¬ 
scriptor for the convergence. Nonetheless, configurations with highly variable 
outcomes usually produce an average that is qualitatively distinct from the rest 
in section 


(see figure 10 


5.2.21. We took the average number of encounters dur¬ 


ing the simulations to compare the context switching model with the context 
permeability model. 
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K-Regular Networks With Number o1 Networks = 4 
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(c) 


K-Reguiar Networks With Number of Networks = 5 



Figure 20: Average number of encounters during a simulation for 100 
independent runs with k-regular topologies in all the network layers. 
We compare the results of context switching with context permeability 
with a numberof networks = {2,3,4, 5}. The switching probability is 
the same in all the networks. 


Figure shows that for 2 networks the results in terms of number of en¬ 
counters for the context switching approximate those for context permeability. 
When we increase the number of networks (n > 3), the configurations with 
higher switching probability lead the results to be closer to what happens in 
context permeability. This is no surprise since more switching makes agents 
switch more often between networks and consequently allows them too choose 
more often from different neighbourhoods. This is almost the same as having 
a larger neighbourhood to choose from, which is what happened in the context 
permeability model. 

The surprise was that for values of A: > 20, the results where practically the 
same for both models independently of the switching probability - even when we 
increased the number of networks, the number of encounters remained around 
2000 on average for all the models. Above a certain level of neighbourhood size 
(and overlapping, which causes the average path length to drop), the switching 
probability becomes less influential. It is still important that we get compara- 
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ble results because we are introducing a temporal component that was never 
explored before in these types of opinion dynamics models: the fact that the 
agents can become active in different networks at different points in time. 




Figure 21: Average number of encounters during a simulation for 100 
independent runs with scale-free topologies in all the network layers. 
We compare the results of context switching with context permeability 
with a numberof networks = {2,3,4, 5}. The switching probability is 
the same in all the networks. 


For scale-free networks, figure shows us similar results. Switching with a 
higher probability displays results similar with those of context permeability. In 
this case, adding more networks makes so that lowering switching probabilities 
displays slightly worse results. This influence comes from the neighbourhoods 
being considerably smaller in scale-free networks - even with values of d = 5 - 
so the impact of low switching was bound to be noticed. 

An increasing number of networks seems to always reduce the number of 
encounters to achieve consensus in both models. This happens due to what 
we called permeability. Having multiple points of dissemination that permeate 
between different social networks enhances the overall convergence to consensus. 
With permeability there are more previously isolated nodes being reached in 
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sparsely connected networks such as scale-free networks. 


6. Conclusions 

In this article, we analysed and discussed our modelling framework for multi- 
relational models of social spaces. These results show that not only contexts 
play an important role in the dissemination of consensus in artificial agent soci¬ 
eties, but also that simple mechanisms can generate a great deal of complexity 
(especially in those scenarios where conventions are being co-learned by mul¬ 
tiple agents). Adding more networks to our opinion dynamics model leads to 
outcomes where consensus was achieved not only more often but also faster. In 
particular, our results show that achieving convergence (i.e. total consensus in a 
network) is not a matter of connectivity. With scale-free networks we managed 
to achieve consensus more often and quicker than with less clustered networks. 
This conclusion is especially important for application (or models) dealing with 
real-world social phenomena, many times represented with such networks. 

The second innovation presented in this article was the fact that our model 
describes an abstract way to represent the time spent on each network using 
the switching probability. This temporal component introduces a new dynamic 
to the study of opinion dynamics. Considering that in agent neighbourhoods, 
agents might not be available at all times - and spend different amounts of 
time in different neighbourhoods-, this is an important factor to be included in 
simulation models (either representing real-world systems with different levels 
of abstraction, or artihcial agent societies for particular applications). 

Future work includes the analysis of different complex network models. 
While k-regular and scale-free networks are two of the most pervasive network 
examples in social simulation, there has been a growing number of models with 
different properties inspired by real-world phenomena that can be considered. 
We will also analyse the networks according to different properties such as be- 
tweeness centrality of certain nodes. Moreover, we plan to apply what we learned 
with our consensus games to cooperation problems. Constructing agents that 
can learn both behaviour and coordination procedures is a difficult task (espe¬ 
cially when working with large distributed artificial societies). Consensus games 
such as this one can be used to reduce the search space in problems where agents 
need to cooperate or coordinate in a decentralised fashion. 

While the convergence and number of encounters seem to be correlated with 
the average path length, this might not be the only property that plays a role in 
faster convergence to consensus. Correlation is just an informal tool that helps 
us assert if some monotonic relationship exists. More research is needed to find 
out, for instance, if self-reinforced structures -like the ones reported by [55]- 
exist in our multi-network structures. In particular, we are interested in finding 
out if permeability mitigates the effects of such phenomenon. 

We have seen that context permeability by context overlapping or switching 
can provide a fairly straightforward modelling methodology. This can be used 
to compose simulation scenarios for more complex social spaces. What remains 
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to be done is to identify more precise contextual strnctures in real-world net¬ 
works. While capturing the structure of real social network can be a daunting 
task, we do have pervasive records of networking activity between social actors: 
online social networks. Analysing these can lead us to insights about complex 
structural properties that emerge from user activity, and contexts that can be 
identified within those networks. A difficult -but not impossible- prospect for 
further research would be to track similar contexts between different real online 
social networks. 


Acknowledgements 

Work supported by the Fundagao para a Ciencia e a Tecnologia under the 
grant SFRH/BD/86034/2012 and by centre grant (to BioISI, Centre Reference: 
UID/MULTI/04046/2013), from FCT/MCTES/PIDDAC, Portugal 


32 



Clustering Coefficient 


A. Overlapping Network Properties 

A.l. Homogeneous Network Topology Configurations 


25 - “ 


20- 



K 

02 

^3 

^5 



2 3 4 5 

Number of Networks 



12 3 15 

Number of Networks 


K 

^ 10 
■ 20 
^ 30 
^ 40 
^ 50 


(a) 


(b) 


Figure A.l: Average path length for 100 instances of overlap¬ 
ping k-regular networks(containing 100 nodes each) with k = 
{1,2,3,4,5} dAlal , and k = {10,20,30,40,50} dAlb} . Since there 
is barely any variation in the property values for each configuration, the 
colours cannot be seen correctly, they are presented in the same order 
as the legend nonetheless. 
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Figure A.2: Clustering coefficient for 100 instances of overlap¬ 
ping k-regular networks (containing 100 nodes each) with k = 
{1,2,3,4,5} (A.2a) and k = {10,20,30,40,50} ( A.2b| ). Note that the 
variation is bigger for networks with less connectivity. Also, for k = 50, 
the network is fully connected. 
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Figure A.3: Average path length and clustering coefficient for 100 in¬ 
stances of overlapping scale-free networks. D is a parameter used to 
construct the network. It dictates how many edges are added by prefer¬ 
ential attachment each time a new node is added to the networks. For 
more details refer to 


A.2. Heterogeneous Network Topology Configurations 



(a) 


(b) 


Figure A.4: Average path length for 100 instances of two overlapping 
networks: 1 scale-free with d = {1,2, 3,4, 5} and 1 k-regular network 
with k = {1,2, 3,4,5} dAl^ , and k = {10,20,30,40,50} ( |X4bt . Since 
there is barely any variation in the property values for each configu¬ 
ration, the colours cannot be seen correctly in some cases, they are 
presented in the same order as the legend nonetheless. 
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Figure A.5: Clustering coefficient for 100 instances of two overlapping 
networks: 1 scale-free with d = {1,2, 3,4, 5} and 1 k-regular network 
with k = {1,2, 3,4,5} dAA^ , and k = {10,20,30,40,50} COb} . Since 
there is barely any variation in the property values for each configu¬ 
ration, the colours cannot be seen correctly in some cases, they are 
presented in the same order as the legend nonetheless. 
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B. Context Permeability Analysis 


Table B.l: Average number of meetings during a simulation run for 100 agents and a maximum of 2000 simulation 
cycles for different configuration of k-regular networks. 


k 

avg. 

1 

sd 

avg. 

Number of Networks 
2 3 

sd avg. sd 

avg. 

4 

sd 

avg. 

5 

sd 

1 

200000 

0.0 

191237 

35780.6 

97770 

86648.0 

44827 

71714.8 

14210 

31633.6 

2 

200000 

0.0 

99744 

89096.9 

22690 

50813.5 

7514 

21352.6 

4602 

6116.7 

3 

198033 

19670.0 

68555 

82808.5 

10217 

33963.1 

6539 

22227.3 

2828 

2283.6 

4 

198093 

19070.0 

38359 

69305.4 

5028 

10680.0 

3710 

4618.7 

2678 

2127.4 

5 

196042 

27845.6 

21287 

48141.1 

3444 

5144.0 

2834 

4334.5 

2278 

1921.5 

10 

178207 

62301.7 

8066 

28298.4 

2894 

6424.5 

1800 

1148.2 

2180 

2078.3 

20 

75541 

95871.3 

2338 

2739.9 

2178 

1588.3 

1936 

989.4 

1802 

823.6 

30 

11703 

43425.0 

1952 

1087.9 

1726 

936.6 

1930 

1096.3 

1956 

1273.5 

40 

1694 

954.5 

1646 

821.9 

1862 

872.3 

1806 

884.6 

1686 

647.1 

50 

1788 

786.9 

1724 

744.7 

1776 

827.7 

1842 

896.8 

1896 

911.4 



Table B.2: Average number of meetings during a simulation run for 100 agents and a maximum of 2000 simulation 
cycles for different configuration of scale-free networks. 


Number of Networks 


1 2 3 4 5 


d 

avg. 

sd 

avg. 

sd 

avg. 

sd 

avg. 

sd 

avg. 

sd 

1 

200000 

0 

141085 

83720 

66020 

82601 

17712 

42444 

14917 

38343.6 

2 

150632 

78710 

21654 

41341 

7555 

21407 

5890 

10889 

3638 

5127.0 

3 

37824 

64073 

12309 

36441 

5299 

20035 

2804 

3791 

2618 

2445.4 

4 

13751 

33295 

4846 

7921 

2852 

2766 

2370 

1441 

2200 

1433.4 

5 

8802 

18754 

2536 

1880 

3204 

11008 

2166 

1269 

1958 

971.6 


Table B.3: Average number of meetings during a simulation run for 100 agents and a maximum of 2000 simulation 
cycles for different configuration of 1 k-regular and 1 scale-free networks. 


d 


1 2 3 4 5 


k 

avg. 

sd 

avg. 

sd 

avg. 

sd 

avg. 

sd 

avg. 

sd 

1 

173312 

61519 

100175 

88406 

63191 

81071 

73715 

84549 

51313 

67704 

2 

118591 

90107 

72054 

86077 

43308 

70318 

31756 

58760 

43281 

72182 

3 

109155 

88722 

38270 

66714 

34899 

65533 

26412 

55589 

23436 

49978 

4 

80035 

88719 

33339 

64075 

21289 

45725 

15575 

39516 

24034 

52407 

5 

89179 

91823 

31044 

62485 

13908 

36791 

14112 

39490 

9252 

28404 

10 

62416 

81036 

11017 

35071 

9095 

27266 

12182 

39798 

5559 

20448 

20 

61371 

81409 

6680 

17900 

6944 

27826 

4721 

19832 

2390 

2313 

30 

51559 

77468 

4980 

12740 

5095 

20028 

2268 

1355 

2046 

1225 

40 

58625 

83049 

5483 

20169 

2726 

3735 

2346 

1503 

2412 

1874 

50 

48702 

74718 

6568 

22611 

2566 

1716 

2362 

1619 

1980 

1081 



C. Context Switching Analysis 



Figure C.l: Perpective plot for the average number of meets during a 
simulation for 100 independent runs: 2 10-regular networks {k = 30). 



Figure C.2: Perpective plot for the average number of meets during a 
simulation for 100 independent runs: 2 30-regular networks (k = 30). 



Figure C.3: Perpective plot for the average number of meets during a 
simulation for 100 independent runs: 2 scale-free networks with d = 1. 
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Figure C.4: Perpective plot for the average number of meets during a 
simulation for 100 independent runs: 2 scale-free networks with d = 5. 



Figure C.5: Perspective plot for the average number of meets during 
a simulation for 100 independent runs: 2 scale-free networks, the first 
with d = 5 and the second with d = 1. 



Figure C.6: Perspective plot for the average number of meets during a 
simulation for 100 independent runs: 1 10-regular network (fc = 10) and 
a scale-free network with d = 1. 
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